(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     70797,       1639]
NotebookOptionsPosition[     69279,       1606]
NotebookOutlinePosition[     69634,       1622]
CellTagsIndexPosition[     69591,       1619]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Theta]", "=", 
    RowBox[{"\[Pi]", "*", 
     RowBox[{"(", 
      RowBox[{"3", "/", "180"}], ")"}]}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Canting", " ", "angle", " ", "of", " ", "the", " ", "applied", " ", 
    "field", " ", "with", " ", "respect", " ", "to", " ", "z"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Epsilon]", "=", "0.01"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"For", " ", "illustrative", " ", "purposes"}], ",", " ", 
    RowBox[{
    "the", " ", "microwave", " ", "power", " ", "characterized", " ", "by", 
     " ", "\[Epsilon]", " ", "is", " ", "taken", " ", "to", " ", "be", " ", 
     "much", " ", "larger", " ", "than", " ", "the", " ", "actual", " ", 
     "value"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]", "=", 
    RowBox[{"2", "\[Pi]", "*", "0.240"}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"Fixed", " ", "driving", " ", "frequency", " ", "240", "GHz"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Phi]", "=", "0"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"\[Phi]", "=", 
     RowBox[{
      RowBox[{"0", " ", "for", " ", "right"}], "-", 
      RowBox[{"handed", " ", "polarization"}]}]}], ";", " ", 
    RowBox[{"\[Phi]", "=", 
     RowBox[{
      RowBox[{"\[Pi]", "/", "2"}], " ", "for", " ", "linear", " ", 
      "polarization"}]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]m", "=", "0.001"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Initial", " ", "spin", " ", "fluctuation", " ", "to", " ", "reduce", " ", 
    "incubation", " ", "time"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tmax", "=", "2000"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Time", " ", "limit", " ", "of", " ", "simulation", " ", "2", "ns"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]E", "=", 
    RowBox[{"2", "\[Pi]", "*", "1.317"}]}], ";", 
   RowBox[{"\[Omega]K", "=", 
    RowBox[{
     RowBox[{"-", "2"}], "\[Pi]", "*", "0.026"}]}], ";", 
   RowBox[{"\[Omega]a", "=", 
    RowBox[{
     RowBox[{"-", "2"}], "\[Pi]", "*", "0.001"}]}], ";", 
   RowBox[{"\[Alpha]", "=", "0.005"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Gilbert", " ", "damping", " ", "is", " ", "amplified", " ", "in", " ", 
     "order", " ", "to", " ", "achieve", " ", "steady"}], "-", 
    RowBox[{
    "state", " ", "oscillation", " ", "quickly", " ", "and", " ", "reduce", 
     " ", "computation", " ", "time"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"h", "[", "t_", "]"}], ":=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"h0", " ", 
        RowBox[{"Sin", "[", "\[Theta]", "]"}]}], "+", 
       RowBox[{"\[Epsilon]", " ", 
        RowBox[{"Cos", "[", 
         RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], ",", 
      RowBox[{"\[Epsilon]", " ", 
       RowBox[{"Sin", "[", 
        RowBox[{
         RowBox[{"\[Omega]", " ", "t"}], "+", "\[Phi]"}], "]"}]}], ",", 
      RowBox[{"h0", " ", 
       RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "}"}]}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"mA", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"mAx", "[", "t", "]"}], ",", 
     RowBox[{"mAy", "[", "t", "]"}], ",", 
     RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"mB", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"mBx", "[", "t", "]"}], ",", 
     RowBox[{"mBy", "[", "t", "]"}], ",", 
     RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LLG", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mA", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"h", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mA", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mB", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mAx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mA", "'"}], "[", "t", "]"}]}]}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"h", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mB", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mA", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mBx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mB", "'"}], "[", "t", "]"}]}]}]}]}], "]"}]}], "}"}], "//",
       "Flatten"}], ")"}], "~", "Join", "~", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"mAx", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mAy", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mAz", "[", "0", "]"}], "\[Equal]", 
       SqrtBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]m", "^", "2"}]}]]}], ",", 
      RowBox[{
       RowBox[{"mBx", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mBy", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mBz", "[", "0", "]"}], "\[Equal]", 
       RowBox[{"-", 
        SqrtBox[
         RowBox[{"1", "-", 
          RowBox[{"\[Delta]m", "^", "2"}]}]]}]}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"variables", "=", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"mA", "[", "t", "]"}], ",", 
      RowBox[{"mB", "[", "t", "]"}]}], "}"}], "//", "Flatten"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"res", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"h0", "=", 
       RowBox[{"2", "\[Pi]", "*", "0.001", "*", 
        RowBox[{"(", 
         RowBox[{"j", "-", "500"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"sol", "=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{"LLG", ",", "variables", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"MaxStepSize", "\[Rule]", "0.001"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Is", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "2", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "1", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}]}]}], ")"}], 
        "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "5", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "4", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "4", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "5", "]"}], "]"}], ",", "t"}], "]"}]}]}], 
         ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{"Is", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"tmax", "/", "2"}], ",", "tmax"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "1000"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"res", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"FrameStyle", "\[Rule]", "Medium"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.784226962268992*^9, 3.784227228761891*^9}, {
   3.7842275027982073`*^9, 3.784227589442974*^9}, {3.784227655152667*^9, 
   3.7842276554218817`*^9}, 3.784227851779511*^9, {3.784228109319489*^9, 
   3.784228151830392*^9}, {3.784228353040759*^9, 3.784228360587185*^9}, {
   3.784228438570054*^9, 3.78422843886664*^9}, {3.784228529783247*^9, 
   3.784228556552698*^9}, {3.784228614972218*^9, 3.784228707442814*^9}, {
   3.784228772223241*^9, 3.7842288038547077`*^9}, 3.7842288833757763`*^9, {
   3.784228964278302*^9, 3.7842290547253933`*^9}, {3.784229105687675*^9, 
   3.784229211003819*^9}, {3.784229301049076*^9, 3.784229538772077*^9}, {
   3.784233261327928*^9, 3.784233290456541*^9}, {3.784233509907217*^9, 
   3.78423351924338*^9}, {3.784233643898909*^9, 3.784233644704205*^9}, {
   3.7842338436015463`*^9, 3.78423384415106*^9}, {3.7842338859399776`*^9, 
   3.784234050734558*^9}, {3.784234117964368*^9, 3.784234220052305*^9}, {
   3.7842346088995457`*^9, 3.784234623687046*^9}, {3.784234664690909*^9, 
   3.784234701419993*^9}, {3.784234898411098*^9, 3.7842349757371283`*^9}, {
   3.7842350472771473`*^9, 3.784235061835257*^9}, {3.784235166316753*^9, 
   3.7842351732323503`*^9}, {3.784235251324011*^9, 3.7842352574040956`*^9}, {
   3.784235289484105*^9, 3.784235327281433*^9}, {3.784235375854424*^9, 
   3.784235376396117*^9}, {3.784235649852745*^9, 3.784235650626583*^9}, {
   3.78423569553428*^9, 3.7842356957866993`*^9}, {3.784236015894725*^9, 
   3.7842360161375713`*^9}, {3.784236053732073*^9, 3.784236101939043*^9}, {
   3.784236156334631*^9, 3.784236207295868*^9}, {3.784236285437698*^9, 
   3.7842362866907177`*^9}, {3.784236327381797*^9, 3.784236385351768*^9}, {
   3.784236520622822*^9, 3.784236541128675*^9}, {3.784236619922337*^9, 
   3.7842366585641336`*^9}, {3.784236759451808*^9, 3.7842368020890417`*^9}, {
   3.7842368533120537`*^9, 3.784236908228676*^9}, {3.784236950202259*^9, 
   3.784237210409375*^9}, {3.784237330445344*^9, 3.7842373643425093`*^9}, {
   3.784242346085524*^9, 3.784242700084516*^9}, 3.784242803084869*^9, {
   3.7842428335804234`*^9, 3.7842429554282722`*^9}, {3.784487858736527*^9, 
   3.7844879036983147`*^9}, {3.784487938035633*^9, 3.784488013832548*^9}, {
   3.784488079260935*^9, 3.784488112647812*^9}, {3.784488232448867*^9, 
   3.784488259390422*^9}, {3.784488307939101*^9, 3.7844884442069387`*^9}, {
   3.784488591758436*^9, 3.784488591946334*^9}, {3.7844886250955133`*^9, 
   3.784488661072905*^9}, 3.784488693129335*^9, {3.784488745200171*^9, 
   3.784488745533964*^9}, 3.784489223120585*^9, {3.78451312066894*^9, 
   3.784513213646172*^9}, {3.7845132715992002`*^9, 3.784513272524568*^9}, 
   3.7845133378267393`*^9, {3.784513472676453*^9, 3.784513477048867*^9}, 
   3.7845769620912733`*^9, {3.784577088884696*^9, 3.7845770969403677`*^9}, {
   3.784577134205518*^9, 3.784577185778541*^9}, {3.7845772229635477`*^9, 
   3.784577261196495*^9}, {3.784577291750091*^9, 3.784577388440825*^9}, 
   3.784577614211357*^9, {3.7846574629201317`*^9, 3.784657463387754*^9}, 
   3.784657500864352*^9, {3.784657532424686*^9, 3.784657540759239*^9}, {
   3.784657594267646*^9, 3.784657617996209*^9}, {3.784657700354312*^9, 
   3.784657744489027*^9}, 3.784657785603743*^9, {3.7846578371437817`*^9, 
   3.7846578377700567`*^9}, {3.784657888400833*^9, 3.784657888817245*^9}, {
   3.7846579255207663`*^9, 3.784657973082341*^9}, {3.784658006523283*^9, 
   3.784658007889015*^9}, 3.784658054151285*^9, {3.784658086438966*^9, 
   3.7846580870994997`*^9}, {3.7846583908921537`*^9, 
   3.7846583918218536`*^9}, {3.7846584357558284`*^9, 
   3.7846584361137543`*^9}, {3.784658741998811*^9, 3.784658849602047*^9}, 
   3.7846588806535096`*^9, {3.784658930815446*^9, 3.7846590110481186`*^9}, {
   3.7846590452214074`*^9, 3.784659048077553*^9}, {3.7846605008772554`*^9, 
   3.7846605089242964`*^9}, 3.7846619574817944`*^9, {3.7846619962482967`*^9, 
   3.784661996888412*^9}, {3.784665048917243*^9, 3.784665059385776*^9}, {
   3.784668036183318*^9, 3.7846680373404293`*^9}, {3.7846680855531187`*^9, 
   3.7846681076236315`*^9}, {3.78467042542705*^9, 3.7846704292615333`*^9}, {
   3.7847012770318804`*^9, 3.7847012985769644`*^9}, 3.7847019007942586`*^9, {
   3.784774322446614*^9, 3.7847743450768642`*^9}, 3.7847775038679357`*^9, {
   3.7848440690584774`*^9, 3.7848440751218033`*^9}, 3.7848446919116893`*^9, 
   3.784848319840334*^9, {3.7911276788422937`*^9, 3.791127940625383*^9}, {
   3.791128486597983*^9, 
   3.791128491939816*^9}},ExpressionUUID->"5981de4b-bc46-466a-bea8-\
49eb3958bb2a"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[1, 0, 0], PointSize[0.004583333333333334], AbsoluteThickness[
     1.6], LineBox[CompressedData["
1:eJw9mnlcTO/7xseehGxJxGTNEpWlFHUVLRSVRIT2fZv2aZ9qppmmZkpIQrKl
KBIR6tMglCzho7IlFeJDypJE+vV7dc/XP17vnrnnPte5nnPO81xnVJz8N7kO
ZjAYH4YwGP///8C/Dv1rwaXTx7qZV9AfcLfFrsfKS8oyqGOf0Ev1k7I8lmY8
m3k7UMoTsbP9+5MfoVJWhMeOlmDlSClPA9sideqqWCkzkXH48iPLBCnPxM70
ORnb+VKejewNzVttk6U8FxNz236ZpkpZFddKsy0W7pHyAlT86zX1b4aUF8FJ
y2CcJEvKi1Fw3+Sxf7aU1eF7cvSXUcekrIHZV9ta0k9KWRP/Lv2pOCRfykux
ZO/TN9sLpLwMW+JGRGaek/Jy5J86PfJqsZRX4LGKSVZ5iZS1UHenaPrRUilr
o0Aj3NfpmpRXYr/VYfSVS1kHDoGzNUIkUtZF6YV7V6/fkPIqKHxdtrS5Usqr
8TxZ79GT21LWw4qRs3L2VElZH5uK2bzpd6UMrP/J0w2oIeYACrNVLAT3pOMG
2Mvxm+5wXzpugCODNz35LmWGIZKbvzobPpCOG2KeWkvKWikz1qAp8sGDrv/V
r8Gsvtvem/9XvxYPc8c52kr7cdZCcrTGtfd/x2eEyhke71ZXS8eNcHLyPz0T
7kjHjTHTaYJXuFQ/xxjfx6kfcb8uHTdBQ2vdiXrp+eSYwOLLHOaNK9JxU1xa
uWzdDKk/HFNM3GjZ+u5/fq6D1+i/uvKnpePrcPHzc4X049Lx9XhR4V7kdEg6
vh7ffh46F7ZXOm6GD5VVyrel85djhv2cxO2G8dJxc1wZKRrfFCYdN8f2mlz1
dB/p+Aaov/R5st5BOr4Bfq+mr/69STq+Edp53yr3rJWOb8TV69zqIcul4xaI
mbL/wMrZ0nELRKmfcpk5XjpuCeX5f+ac/ms2wLDE/twFomsfiDmWUFyq/F7/
CbHEEg27AgfPuEbMsMKoTaMm6R+V1lvBfvDRLC5PWm8FvU31yx+7S+ut8KvL
UChjKq3fhI/asZOGz5XWb0Lpty9DTg2S1m/CEfvgIbeer6f6TYjtGjFK6zwx
wxpZb70XPOISwxrCCedMXWyIOf0cY99cPlNab43GYxmSU5/WUf1m/Dnhfvp9
MTE2gy/qktUPIeZsxrINp/d4LSWWbEbeKt/HKp9Nqd4G4rVxVxWOEcMGN5z3
qv6wIubYwGbbamunPyZUb4Pbl+KX/MkhZmxBpefg27b6xNgCn1t7AobVG1P9
FvTpdI1JdCeWbMGNIKH85g4jqt8KJT/riFcsYmyFeex4SULbWqrfCot/E+Kz
thBLtqLqUtJct6trqN4Wm/QHGypNIIYtSgVnGurtDaneFonmO53OHDGgeltU
cdftFj4C1W/DMv7V/ONt+lS/DSnjskcljNaj+m0I0toZ+GPMKqrv/zz8Nv78
sJLqt+OujuJCrYNaVL8d33VbPORmLKf67Vgve6PQJ0CT6rcjjz22tGLqEqq3
w4mG2WmtjapUb4fFt8perFOZSfV2OBRkv9K3cBLV28HBOFTXI+XnPwP1O2Cq
UdybETJcf6B+BxbJOUxZ80B5gDk7sEBgY3jxr+oAS3bA50LR0Rt/NAaYsRN+
Wcl/635qUf1OLHYdWaYeo0f1O3Hp657UAw8NqX4nVj9qZM+/ZEz1uyBA2YeS
P+upfheMr17ISz9rQfW7YHh74YftPpuofhe0U+7PYxjZUL092scN9bjraUv1
9tBYdo+X02pH9fZ4ydR981diT/X2UK02K5g53InqHfDpITthlLfLADMdUH9S
Y5GHojt9nwPWd3dXyE3xGmAHB2xMHl34INyXvt8BJpMcLv4+wRrgHAf8GTXu
9sf7gdTPAX/j/rl81TBkgJscUHAouO22PJv6O4LR5l1uaxpJ/R1xYl9Pql9v
DPV3hOIPm10OJ+KovyPeBCimiWZzqb8jWjoqVfaN51N/R7T3Jip9TxRSf0d8
n895+y1HTP0dIV/hcDmkcjf1d0LugurQwit7qb8T9kZN4B3elEn9nbDlakvp
n7JD1N8JCkoteR1KR6m/E7oP15Vc5J+g/k4YHNNqOOVlHvV3gvZ5j3s6Tmep
vxN8umyOV5YXU39nKIXfUpxndZn6OyP1p/2kKerl1N8Z23lTZb+su0H9ncHz
OLfbsew29XfGeNlM52c5NdTfGSmN7OeOE2upvzPc9XPzU4WPqL8zttVMP7Xo
LDHDBYxXxjX/cOjzTBd8fcfSODn6PvV3QWjobrWLxVXU3wWpy9b+J3+ukvq7
YMSZ6vNl665Tfxecnlnw0nY0Hb/EBUb10Z81ckqpf3+9j8e8qQol1N8VfV1d
6bZ6dD6YrvB/t9vlVPM56u+K7fqrTcTMAurvikCno+NdRuRTf1fMunc6wyso
l/r3f17+pJLrOvJD4oqNvwIKg4KPUX9X6N7p8AlrzaH+bvB3e5Bwa8oR6u+G
B/F5e4vvkN9wQ/77zjuTbmZRfzc4WQ9TCuqj+cFxw9fvCzPKAvdTfzc4rNHx
qJ6YQf3d8FokiW9vpPnV5IbCrgvCczV7qL87zB0VreLr0qm/O2ytDYuP/Kb5
CXcMthSMeq1F7OCOCbs2WB0QpFF/dzA+Hjk+uz6V+rujcn9f6P1ymu8Sd9w+
q6r5+7qI+rvjUOVc0a+mFOrvgW7ReqtYRWKmB2Tn2/sdcUum/h7weKXQV1ZN
15ODB7oMeG1xIOZ44MGKV5r11UnU3wNF81cZv3YhlnjgHjumtVOeuMkDqyy2
zD5wX0D9PaGwre9VahYx0xNNd5Mdr4cRwxOnB+83WO9K7OCJHzXtzyOciTme
mHjF7fjzQOIcT+y4NblHbS+xxBPG6dW86CriJk9Ehab9nDGOjofhhUWZCbZ3
/ImZXlgX5vhs9ztieKFnEbxnsKX6vZB7XKPh7Dw6PxwvdB+3um74izjHC/OC
z9o+7qTzKfGCrPBwqnCy9P7jhcyD1VNnB5JfDG+oV24pPcsjP5neKHq5XOdu
mNR/b8xMFLRt30/zw8Ebl7U8N1xup/nD8YaJbePlmOR91N8b3fGPO5/F0nyU
eOOX7I9/X7ceoP7e2JCs3Bj9kuY3wwejfS81Tzan64HpgythQ/98PkfXC3yQ
9bIkvO3BServg84Rn+IX9EqvPx9crgwsNr9O12uODxStjGQPD71I/X0wdwps
orWuUn8f4JrH4tYRdD9j+CKWZb9zSOxd6u+L8krcCl7wlPr74t8XWtuenmqm
/r44ts7A/fH8TurvCzPV9uSdt39Rf1+cKq6YNkirnfr7ImD+juJAg2fU3xda
fWt2jZlcTf39ICnam5e/qoz6+8FH3nmOzfwz1N8Ptn4NHR8PHaf+fvhqFTX8
99CD1N8P2+7fj1Z/s69ioL8fNBU3bPwiGzTAEj9Yp7e5X+6l52OTH1qr7+gM
lhcOjDP8YZly/u1yRsoAM/3hmjHY45nc0QGGPzzip0S9iLowwA79n+crpIlY
xwaY44+Oawi3kium/v7QmP8oPTf7BPX3x4ctYY0Jj4oGuMkfciVuBc/1z1F/
FoZyni9smpg/wPIsmAv3HXdUK6DjYWHNylYbW02qV2chRn6uosxwYrDwtPbo
DLbZ2QG2ZGHFrDz/5pJCOl4WGrKO3pokpHEWC9fMGriH26g/h4Wfs8w722Pp
+9JYMO5YcMG4hjiHhex9AT1nPhEXsfBu6dkf0wefJ30spLtdeaIwnbiWhdOu
PzdfMiZuYsFdfmVzYgRxR3+/KXYjR5QQMwKQ9Pe1RWEnsXwAVrpyHK3V6Hwy
A8B91zM4yI1YPQBOg5yOVB4iRgA2qvAvlNQSWwYgNvx1/s/BUr8C8MyWOVJn
GTErACnPX0PWlZgTgLF29u9r9hGnBUD8vDuh/DZxTgAOPVZrG9RNXBQAVuHZ
4vz5F0l/AGxVnoy/ZkdcGwDrM8cstMTETf3jfW6m068TdwTg/aOprqzvxIxA
iB07gpaolpD+QASeOP1q1w5iZiA+lU+PGbKbWD0QvUuO28y6TYxAuNRZad76
TWwZCNGFXscOjUukPxCVzvVBBR7ErEB4ynXIDsoh5gRCedZYyY964rRA5CW9
G5suf5n0B+LC63ObXq0jLgqEx1DNq40JxJJAzFo4b/Phf4hrA6FnajSW2UPc
FIiCHUHXfVeUkv5AdG/IXSwKImYEgXX1yFjeeWL5IIT6sRR3dBAzg8C5u/uv
ovoV0h+Et9NNYm6wiBEErTkSe/tiYssglL0VR3Z9J3YIQnHtw3/E2ldJfxC0
5fsWL4km5gThxZygz403iNOCsC6iWOfkyGukPwjKrNTwJCvioiDwqhdp7ski
lgRhzNLWiXdbiWuD4LpnXu9y9TLSH4TA9ebfG6OIO4Jgs9xj2f27xIxgjL8k
u3CkUjnpD4byfqWog17EzGDo3G8YKiojVg/GNu+Tje/G/kP6g7HUZKhTsQux
ZTB+RV+a/+sqsUMwynWeT7o1voL0B+PYUGULpg8xJxi7upjrJ90hTguGz0lv
1/pxEtIfjMpjE77/BnFRMGR+Vo74yyKWBON+QjC35ShxbTCsJptGFz8hbgqG
2K5wRcDw66Q/GAuvR6rO1SFmhKBltBajwY9YPgS6DTuOi04QM0OgJWrS2PiC
WD0EW9MuyqhMuEH6Q/BvfFjlaHNiyxCc6NBaOy2R2KH/+yZ6Fmy+TswKwT6X
L54VvcScEMzbZrPcVfcm6Q/BDl81BYsI4pwQZEky3BKvEheF4EpZT93EP8SS
ELQXfwwYql9J+kPgccswxz6BuKn/+7ax1DSqiTtCoH9vuDlX/hbpD8XoIbFZ
LtuI5UOxRtVItfE4MTMUybaDI3u/EKuHYuhunb0PV98m/aFwCVrw01FEbBmK
W+Nlz91oJHbo5zcHbRiad0h/KH6cVlFawCfmhKKzqaDKopE4LRTNzz88itOq
Iv2hcNC0tn6UTlwUis/T3syw7CCWhKJyuX61vGU16Q/FBUVR/txi4qZQpCSd
+3xS4S7pD0VPwuac49HEjDAwBnNLlr8nlg/DpbJhk4Osa0h/GGYeCT8fdINY
vf/zdyziLZfdI/1hUBm+N3dWPrFlGB5fvOYlo3Kf9IchLDLIbuohYlYYfEzZ
WyOUHpD+MLDfFB9be4g4LQyvJKk9WSoPSX8YDumeO3PkNHFRGHY18Ddfm1pL
+sNwz9r38UwX4tow/G7xmpR1lrgpDBOX2DUu/03cEYbKZau6R61/RPrZiD/m
Y7r0ILE8G8wLCk0324mZbGiEPh9Tv/Yx6WdDM/dlKP8wMdg4sGPbzM5uYks2
JJ43sjS2PiH9bJxXuTjIo5SYxcYa/8qFJdP+Jf1svIoKK9flEqexYX/Ydz3z
C3EOG42vh4Ul7XxK+tm4fOTutJyHxBI2Fo2pCEg1qiP9bHiHdmjyKoib2Pjz
X8rj86vrST8bp3a05ZtKiBnh2MY6v4Rj0kD6wxH6pkUz/QkxMxybXjYolbk8
I/3hOJDwd5H+b2KEw0v9RYPD/uekPxzOAW2y3itfkP5wqJ+U4Re8IWaFY2th
tFl42kvSH45PKsq9C4xfkf5wVOp17lQZ0kj6w3Gn1ejhjTvEReGwjNv5Omrv
a9IfjjFuF0OalzWR/nAIiw46yZ0gbgqH1lC/a32T35D+cJSZD5MdtpuYEYHv
qeFWXmObSX8Eftwflpqwh5gZgZnjf5ieUW4h/REw83NUW3OWGBG4fObyjYtG
raQ/Am1bAhYFthI7RCDguGD4T+Fb0h+BXJPH26tXviP9EcjaWHohopM4LQKf
HnlY5he/J/0RWPxJN3ZQTBvpj0DXz7JZ+ps/kP4IxB6667dkxUfSH4HgdJXa
4tn/kf4IbH6ts6CC+Yn0R0D0j9smp4WfSX8kvruHfKg0aif9keDqPsxSC/hC
+iNRqenjGb+1g/RHQs5Dj+PD6CT9kbj24UyZxkViy0hIHDzG3An6SvojcU7z
6X+r8I30R0LZ0U7xptJ30h8JR5wbUTj0B+mPhLx25O/QPuKc/u+vcu6rkftJ
+iNh8aKtPWVxN+mPxAqr21cb3H6R/ki0+msvUSjtIf2R2LrOOaVv9h/SHwm7
5rrNRUW9pD8KCnuNHpbs6iP9UfgRzPHefYAhGdAfBWcZ8WvrA4MGWD0KpkL7
I1vPDx5gRKGlMuuhpHXIAFtGYWIl/6WW1rABdogCb5NO45DC4QPMisKJF4qq
C81kBpgThQmcmlMHJ8sOcFoUtmbFaO6cIjfAOVEoTZ/CqbUdM8BFUfinJ21p
n7/8AEuiwPRZ9jf61bgBro1Crr7+/WVREwa4KQq2J+tNZW0mDXBHFJTKzYqa
AycPMCMaL3WGuc5qmTLA8tFYbue6bvulaaQ/Ghcyph0QMJikPxryrrfTvdpV
SH80pr3R2Hd12WzSH419lbrrtN/PJf3ReDL352+t6QtIfzQyuF3HNX+qkf5o
PFvpnWW3XoP0R8NIp6bOfcwy0h+N+RPOWslu1yL90XihdFl++Upd0h8NU4er
Jhm5+qQ/Gmu3/rn/fJUh6Y/GQ/m6V2p5a0l/NHqK9RvUJMakPwa6qbyOBwqm
pD8G+tV840cbiJkxGGVovGz4BhPSHwPRAq+6chsj0h8DeVXv7/wX1M8yBmdM
Ooxr1oH09/M8p2/Gd+l4WTE4aJ486Egm6eHE4Or4refHjyK9aTG4VjTd4Owk
Oh85MXBznWT725jOV1H/9/klxF5ZNZ/0xwBPRu/its4h/TFgf65cZaU8i/TH
YPkk90tTqsm/jhi8WvqP0bI8ZdIfixVPZZ6XblUi/bFYY9WWfZpL84MZi/lj
NQtOK9L8UY/F7N9Tj9QPovmFWNhnLzN7qkHzzzIWQw8yXk/SGUv6Y1ErjF47
JZjmMysWXzNlp20YTPOdE4tkh5xdExtGkP5Y+NS+LS77SNdPTiy27ZsiuqAx
lPTHQun8mF1JJ+n6k8QionDWkV59uj5r+/vtPyvw/knXc1Ms9Ldsq56eR9d7
RywOXy+060r6Tdc/B9tC7jRX7KP7hwwHcvP0JKNv0f1GngOTwS/K5o/vGmBF
Dua90WnVYNP9i8kBttXbxX6j+50qB7GRFts84ul+qM5BdkPPPf+5dP/U5mCV
XvnBrLt0/wUHMft3K8vupvu1KQfXjiT+9+0C3d8tOfiqMsFEfRXd/205MAvX
U2uZQc8HBw7ej39iO96Ynh8eHDjtHHe2N136vOFAt/im6bhOeh6xOfizUbbg
njUxh4O8RRvNHpbQ80zAQb6Jxfr/JhGncTD5rXH+u2B6HmZyEDKl8+Kif+n5
mdPPr79lL1pOnPf/+k31lQ7Q87eIA/tum1uNDOJSDqZtd/HRfU7PcwkHiy8x
POLK6Plf1d8vrlVDJp/WB7Uc8PoO6w3NofVDAwc384Tfwo7SeqOJg8Is7nOD
QlqPtHEQ7b8vcN5NWr90cLDIqHX9q2Za73RzsH6E/LDRctL1Txx2C2IrtunR
ekkmDu8SX5o7h9P6Sj4O69oPJlaV0XpMMQ4KhW++q8oQM+OQe6X75xw7Wt+p
xsG5cdJe5xLpejEOZ/YN7xYoEGvHYU9qQy1iaH2JOMTeW8te+YnWp6Zx2NL5
ecQcB2LLOJxIkSsOvEDrX9s4hDLlPtmfp/WyQxx+T1k/LOUKra894nDDbK8O
u4bW46w4yB3ewS54R+t3dhyK+CPlL8sSc+Jwz6LUa5EWrf8FcXjmE3kt15v2
D2lxSDVZWHzrFO03MuPwK//bksX/0X4lJw6+i9b+TVpOnBeHr5e9BPxE2v8U
xaHbcHXms1e0fyqNQ+cdGQUtXWJJHJ5E5NRsz6b9WVUcji+oWTBkhHQ/F4el
pbcVvgbTfq8hDnaDhp3/9J72j01xOHCtXPWaPXFbHHoepIhUX0r3u3Go8Mg7
MXgncXccvKpvS6Y00/6YEY9AC59D8CaWiYf1o2m+mR20H5ePx8l3D5zmSffv
ivFI0ilya9Kh/T0zHnOsvfaES/MC1XhM6LqSVZBD+YJ6PIqvTE9f8YryCe14
LJ2bfqN+BjHi8dvGKNLag/IP037evuK9/yVpnhKP+VeXv/0xktg2HjXrBM1F
TpTPOMQjJejoZT8J5T0e/f2+z+X3zSRmxWPFoars5ULKl9jxcMgft/3tD8qr
OPFol3s//K0rsSAe/03Oq5F7TvlYWjzSTR6LV2wizoyHPKY0rX4gzefiUT9c
0ty3gTgvHkOU16+0e0R5YFE8jBimX5fYEpfG40vKklD7Zml+GQ+1cTtC7/kT
V/XrT3r7cjtDmmfGY3eBe/mXPZR/NsQjJ/ukccACaZ4bj3WG7Z4iPuWpbfGo
8/1XWFRBeWtHPOSazayn91Ie2x0Ps4v13af1iRkJMPQO9V3Ep7xXJgFzpo01
i3hyhvxPgPb5tgrf2cSKCSgyMo56E36a/E/Ar+znOeefUH6smoDG8y6F1zSI
1RPg6bYquHFPHvmfAI3wmOE9v06R/wl4LJmR/8WZ2DQBywUlt48+yiX/E5Af
d6L5myGxbQJcdEeI7l4+Sf4nILJzztZR6sQeCTgiCq86dIbyb1YCWjtNNrst
IGYnYGH5JQvLguPkfwL2ZW9tW6NBLEhAkN+KqQuvUL6eloBNRw0LOtYQZyag
pWDJdM4jyudzEqC0qPj0LUfivP7zaaSge/Z7DvmfABtTXcd5QuLSBFjn1VrO
mEksSYDtlffnR4YfIf8TkOswrUxcnE3+JyDm1J9ThR2Hyf8EdHJGB5gsJW5K
gFfyp6l6EYfI/wTcW61VEHzrIPnfX//6itLDicTdCfi8Z6W5pkcW+c/FLf1E
GW7FAfKfC58fcl9PKxHLc/HibG60ODyT/Odi0OE3HxVe7if/uVCcLlylZkis
ysVwr/0a185kkP9cGFWb6Z1RJNbmoj3yJq9NQO9LwEVm84Enu/7sJf+52CMe
ldETSGzJxTy/uQ15n/aQ/1w8zbdetdWT2IGLvrFDf35rSyf/uRgy9mq1rzcx
iwv1N4p/L3TsJv+5aFTe6HiJTczhIm30QX3XIcQCLqy/DHY4n5ZG/nPBH7WU
kcIkzuTCqlvj6H3XVPKfi9HLp75xTxWT/1ws+eL1tKZcRP5zccZLxvl1B73f
KeUitPe6VpIqsYSLAN97N4tcksl/LhjMg0v1TtL7oVouZrU2eSp9TCL/uXi0
9NWvxUuJm7goKKlSceMIyH8uFHYn3zr9iE/+9+ube8azbS5xNxfBVeaLR3ES
yX8eWGNi7/a95JH/PCwIOMu/sJpYnoeQ55qrJx7jkv88uNopmMvLEjN56K3o
ycoISSD/eZDNmJ2e0RpP/vMwcs7RyYytxNo8rBp8IPnmvTjyn4crGyWVj4yI
TXmY/uqZ3LQbHPKfh9ByR+W9BsS2PEgWntS27V+zDPjPg8YOpmbhnhjyn4ch
lXXjR5+MJv95aH2WJdpVFkX+8/DBiDk56Vkk+c+D6OSSN8F/Ish/Hp4vNqqT
nUOcxoOK3K99y6zDyX8e2pW2BL5MZJP/PPhJ9Pw+V4SR/zx82ZK42vpvKPnP
Q6Fjo3GPAXEpD8qLl+c9FIaQ/zyMWKCpd70hmPznYVml8rnyhcS1PFw0mh5x
MYHeJzbwEGtTLpvRFEj+8+DCWfjI2pC4jYdBIcq7mvICyH8ehIs/cJZMJO7m
ofP5M+cVXBb5nwidFjXD+k5/8j8RVVduJotb/cj/RFT+4nnxmnzJ/0TYRQec
LnrrQ/4n4kzJ1sfDvnuT/4kw7j/F7FHE6om4X237rXuhF/mfiPjoC7e9bDzJ
/0RY7LnYcpTvQf4noibp5/ytN9zJ/0RsPbknX3kYsW0irg8dmbrFwo38T4Sw
WNZ5bo4r+Z+IA4VNxf5dLuR/It7YNl57ZknMToTYJO/Bp0Jn8j8R6J4NKzli
QSISB4+4oO/tRP4ngmnn176kxpH8T8Qu/11z7i4mzkkEVza5a/ghB/I/Ecp7
J4brKhMX9fc3n1X6N96e/E+EduaJW3UriCWJ+Kcy/0p4OHFVInY/8Ztpq0b1
tYmw1145eekr4oZEGJ4uXfJiD/Vv6vfv5tQn798StyViLVo2dOgRdyTi6+K4
iUYzd5H/iTDo1fyR6bOF/OeDNXjR2TF/6PeNMny0zN3R966Pfp8nz8fHOj/u
3wT6/aAiH8cDD+SbrdhJ/vPxwfZUXkme1H8+GIpRZ4ZPo/mkzofiAts/RXpb
yX8+xJq5sWs+euoP+M/H9ge/xUlfggbYlA+DzxFKXRI6fks+1nLvnGp8Qdej
LR/aZW8yDoijBj7vwAc7U9VfdzL9PsODj4q2K5Ozh9L9jsVHKh5lTVpBv79h
8zHqVNaayqd0/+Tw0fn0fK1aAf1eQsDHAtnti2WX0fMsjY+eXyFyG2ZJn/98
uE/5/c+JDFov5PAh++XQsY96tF7L42OTfon2h0/S95V8TM/Rc1W+T+vbUj4W
7nwnLFai/YaED96NgjlWy2j/VMXHnrX2webrpHkuHzcimK8rjWi/3MDH+tZj
Fdp50v09H3Plx2w26qT8o42PIDtJ873fUwe4gw/NLa93bvxIeVl3//m9/Prj
lbm032cIYP5n2J0jRylvlBHgRcipfX3vieUF6D52JtZzGuWBigIodh7rsjGj
/kwBHho1/WkaTvmmqgB5vK6EITMpP1UXYEn2YV+dsZTPagsw/1Xc+zlbaL8O
AY6z769S4tJ+21QALY2W1jWbaT9tKcD+61qLDmpS3m0rwCe1M5HrvWh/7CBA
9ZzGOtc5tB/2EKAkXDZ9ujHtf1kCOAU7+766RPtbtgBMc7neU+60v+UIMM9Y
7QnHnPazAgEuubnUGewgP9IE8Dj8yKkhlfavmQK4eB4Rlb2h/WuOAF3jTtmm
bSDOE+Bj6sETzx7S/rVIgBnBWrIPXYlLBYjW2huoOVr6/kaAOtaLKZ4C2t9W
CWDr1W5XMZ64tt+PTR/Mx16m/W5D//k+XJYUEUTc1H+8600MhcbEbQKYVtus
Yi4h7hBAXa9+vs8i4m4B2pMZgS2riRlJeJAUulPXmVgmCZlqRrqu0vdR8kmo
Kzn/X/IHYsUkOE/6VjLBjI6PmQQ3N2Vx1nVi1STMcHIy6Z1L+tSTYJIUUJvL
JtZOwkGLg67hdcRIwqsNBwrWG9D5M03CaeVDrUpXiS2T8L3zYb2SIZ1/2yS0
PXH8Vvyc2CEJV8af8krhkV8eSdjZlaV9aC35yUpCiX/uIg0l8pudhMRWxsRh
I2k+cJIw50BIzaaJlIcIkuDy6OqWTn2aP2lJODV976vQNMpXMpOQ5jKux2UU
5TE5SUhdYG/pdJnmY14SMqabyvMOSd/XJOHJA6fy7U7S/CcJR3YqTFbMpvkt
SYJ5hoHr1JE0/6v6+2dmrmQUUj5V26/v8ZCvQzLoemlIgs7DfVGdTygPa0qC
sdmFjOBIys/akuChX58tl0fvMzqS0P6F1xetT3lcdxKayvllHw/Q9coQQknG
amQTn/I9GSFSTPIVb+rR+wN5IWb8aAlTG0vXu6IQPRaKC6quUf7I7B/3WM80
mkX5paoQFnZVv+57jx5gdSFU9r06tH8z5aXaQoi9cooU6qYMMIS4/fdna9dy
ymdNhTj91ODX2LOUr1sK0bjlg2evObGtEBUbjk448e+MAXYQ4nD6iYL7WpTX
ewixop3r/NFLcYBZQkxImL1HpXriALOFmK9s5WEWR3ktR4ifsbK9jzbQ8QqE
kPWtWfCLTe8n0oR47O+dXh9M+WumENtLbodpqFLemiPE67z2kLWf6f1KnhAF
dg979yfR+5mi/v7Z3YdavlNeWirE6m/rzHYMp7xUIsSqofsvcFbS/bKq/3ie
Llr1Xof8rBXCmRMypuYd5Z0NQigeN/zUPZXyzSYhbsQ5KptLKJ9sE0L/toHp
7Aqabx1CdOjwVGzlaT52C/G06Njdqmk0fxnJ2GCa9eSBmOa7TDJS5rJuFNrQ
9SGfjGeD3u1vcJbmf8mIK1ziZ3SWrj9mMphz5d3aFtL1qpoMbvSJ65nS99Xq
yVitrqrEM6V8TjsZHjculhvPo/wNyUjVzn2zUo/yNtNkGAyP2T45kfI0y2SM
zHjQW9JNeZltMuot36Y3pFE+5pCM9Rf0dg1bT/mYRzJ2lZvHuzKlv7dIhin/
hedpecq/2Mm40Gkl3z2Z8i5OMmYsu7g2153yLEEyRslYZM6upPwqLRmM+8bD
TxhQXpWZjEs7ws78+5LyqJxkDP5y7dm9dMqj8pKxyke+scdB+vxPxrnkg38U
jSh/Kk2GdZWrxqzVlDdJ+o93T0vuRhPKl6qSsWTUS8vBTpQn1SYj2yr8VG0y
5UcNybArNV8TeFP6+7ZkaIcXr7k+grgtGaxiOeOSrZQPdSRj+r4nXytcKB/q
7tf/03/Vv2qUDzFSUPZt5HOVkZQHyaRg5qCGfSY/pPlPCpQiQy8N7aS8RzEF
i+8aeS3/TfkOMwVmjUPuPhwvzX9SYHvZVfnNCsp71FPgEPxNmOxG+Y52CvrO
sLqGHKU8BynQXKU4P+0t5TemKdj8bO0i3lJiyxSEnjN/YiSk/MY2BT7HGX4G
H2j95ZCCwi0GR75aEHukIOrhIuvGclqvsVLwML++Nl2TmJ2CSWNLWnILKa/h
pGDprDiVK4uJBSmoHft6mX8J5TNpKdj4kmEaYkCcmQJoPf3cc4LympwU6Ojq
HY2ppbwmLwVF+45Utw8lLkpBvaHW5fsGlNeUpuDgEOMluYmU10hScHxu6Q3Z
J5TPVKVAnFs3f99c4toUcJpe3XzLobymIQUKzfFyx99QPtOUggdHDqmGmxK3
pcDt46CIKSWU13T0n886g9U684i7UzAutWtxRjatbxkiyHFMV35SIpYRIaP+
w7y+g5TXyIsQk7bWaZ8KsaIIUycuPsM+Q/kNU4QWs39aQ1YSq4pgywh+vaOG
8ht1Ee4ctp3BcCTWFuHqw7wLy3oov4EIyxlhso8ziE1FmG29ofr8CmJLEezi
DuSde0Z5jq0IW8pSrfJjiR1ESLCZox+jKs1/+o+HXX907L+U57BEqFsxD4iX
5j8itKbYu3zXlOY/IjgPvyAZ947yHYEIwuvntscelOY/Ihxc0Bk6x1qa/4iA
x1sSP8sR54jQq2k0n7ub8qA8Ed79/umrPJq4SAS3M9f6pxnlQ6X9x4fbLxLH
EktEqN+odXVMBuVFVSLIDlbp/TGDuFYEhxHme6cVUl7UIMKNjZM/Oq8mbhLB
L3NL1dlayo/aRNhh/t/GN27EHSKYRZm8aOqjPKlbhK/PFH/tOSj9/bEYj4eN
tGpeSSwjhl529pGLz2l/JC/GZker7O/RxIpitHdn/EmZRcwUw/BvtpVnDeVP
qmIMEVTJh4YQq4uhc6kjNkeFWFuMzwFT8uoe0v4MYhj/17NpBIfYVAyl7Fu/
VTSJLcUIY827p/CO8ipbMWovbhXVHSR2EOPd7U8nN1kTe4jhoPZ7VLQcMUuM
hNKnpqZ3KM9iixG4/5v7qQRijhiPLjX9EBkQC8SQdO749YVBnCbG7Hdv2q9f
p7wrU4y7cRYV3QnEOWJsm9lTkWRCnCdGx2t2oZcccZEY6oufH9/9mPKx0v76
4Tum/DhALBGj0n9Hc7ITcZUYbuN5czYsIq4VY0fzum3qPyk/axBD7Ht8jdpN
4iYxhqWe7NJNI24TY1CQ2t3Nu4g7xFg63tWWpUbcLUZe5sZ4fi/lbYxUtJQ9
dtzzgFgmFab+U6xFOcTyqbBRGffeL4hYMRUvZFQ/aJgQM1MxcU6hbe1UYtVU
7JMr91rbSXmdeioi8pfECe4Qa6eibvpczf3ZxP2787JVK1Z5hxKbpqJ792dR
30ZOxf8BUR8stg==
      "]]}}, {}, {}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameStyle->Medium,
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 1001.}, {-48.60794908170873, 1.6109290567424235`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.784851362636764*^9},
 CellLabel->"Out[11]=",ExpressionUUID->"82d99a53-3f41-4f55-9d31-2bd4476e5722"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"res", "[", 
  RowBox[{"[", "481", "]"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"res", "[", 
  RowBox[{"[", "764", "]"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"res", "[", 
  RowBox[{"[", "841", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.7846669847734084`*^9, 3.7846669850828342`*^9}, {
  3.7846670374103565`*^9, 3.7846670403123174`*^9}, {3.7846670771053624`*^9, 
  3.7846670995257454`*^9}, {3.784667144872143*^9, 3.7846671645255747`*^9}, {
  3.784675065126302*^9, 3.7846750701958647`*^9}, {3.7846751613624377`*^9, 
  3.784675183889835*^9}, {3.784675544305384*^9, 3.784675585891162*^9}, {
  3.784692253956747*^9, 3.784692285200254*^9}, {3.7847063523919415`*^9, 
  3.784706460390225*^9}, {3.78477725656559*^9, 3.784777384235519*^9}, {
  3.784780404335057*^9, 3.784780452759684*^9}, {3.784780492787293*^9, 
  3.7847805130770936`*^9}, {3.7848476234870014`*^9, 3.7848476236318283`*^9}, {
  3.784851539518652*^9, 3.784851546284134*^9}},
 CellLabel->"In[15]:=",ExpressionUUID->"b6431c20-3ae4-4860-ab7e-b875fe6ca575"],

Cell[BoxData[
 RowBox[{"-", "48.60794908170873`"}]], "Output",
 CellChangeTimes->{
  3.7846669870806727`*^9, {3.784667089778704*^9, 3.784667099710784*^9}, {
   3.784667147767415*^9, 3.7846671646949863`*^9}, {3.7846750672644787`*^9, 
   3.784675070459072*^9}, {3.784675158738233*^9, 3.7846751839901667`*^9}, {
   3.7846755566311936`*^9, 3.7846755859915*^9}, {3.784692250096774*^9, 
   3.7846922854230986`*^9}, {3.7847063609764595`*^9, 
   3.7847064605308475`*^9}, {3.7847772529194255`*^9, 
   3.7847773844386225`*^9}, {3.7847804006169972`*^9, 
   3.7847804528533955`*^9}, {3.784780493584159*^9, 3.7847805131708302`*^9}, 
   3.7848515474165254`*^9},
 CellLabel->"Out[15]=",ExpressionUUID->"ed1c2fcc-0ad6-41e5-83e7-f66d6fb17979"],

Cell[BoxData[
 RowBox[{"-", "10.971232058254895`"}]], "Output",
 CellChangeTimes->{
  3.7846669870806727`*^9, {3.784667089778704*^9, 3.784667099710784*^9}, {
   3.784667147767415*^9, 3.7846671646949863`*^9}, {3.7846750672644787`*^9, 
   3.784675070459072*^9}, {3.784675158738233*^9, 3.7846751839901667`*^9}, {
   3.7846755566311936`*^9, 3.7846755859915*^9}, {3.784692250096774*^9, 
   3.7846922854230986`*^9}, {3.7847063609764595`*^9, 
   3.7847064605308475`*^9}, {3.7847772529194255`*^9, 
   3.7847773844386225`*^9}, {3.7847804006169972`*^9, 
   3.7847804528533955`*^9}, {3.784780493584159*^9, 3.7847805131708302`*^9}, 
   3.784851547432148*^9},
 CellLabel->"Out[16]=",ExpressionUUID->"036d09fd-b7ef-4dd2-9654-e832b4773e44"],

Cell[BoxData[
 RowBox[{"-", "8.839965619360402`"}]], "Output",
 CellChangeTimes->{
  3.7846669870806727`*^9, {3.784667089778704*^9, 3.784667099710784*^9}, {
   3.784667147767415*^9, 3.7846671646949863`*^9}, {3.7846750672644787`*^9, 
   3.784675070459072*^9}, {3.784675158738233*^9, 3.7846751839901667`*^9}, {
   3.7846755566311936`*^9, 3.7846755859915*^9}, {3.784692250096774*^9, 
   3.7846922854230986`*^9}, {3.7847063609764595`*^9, 
   3.7847064605308475`*^9}, {3.7847772529194255`*^9, 
   3.7847773844386225`*^9}, {3.7847804006169972`*^9, 
   3.7847804528533955`*^9}, {3.784780493584159*^9, 3.7847805131708302`*^9}, 
   3.784851547432148*^9},
 CellLabel->"Out[17]=",ExpressionUUID->"cb8b5f97-1376-45f5-ab4a-9d35c5b21037"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Clear", "[", 
  RowBox[{"h0", ",", "\[Phi]"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h1", "=", 
   RowBox[{"2", "\[Pi]", "*", "0.001", "*", 
    RowBox[{"(", 
     RowBox[{"481", "-", "500"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h2", "=", 
   RowBox[{"2", "\[Pi]", "*", "0.001", "*", 
    RowBox[{"(", 
     RowBox[{"764", "-", "500"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h3", "=", 
   RowBox[{"2", "\[Pi]", "*", "0.001", "*", 
    RowBox[{"(", 
     RowBox[{"841", "-", "500"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"H1", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"h1", " ", 
       RowBox[{"Sin", "[", "\[Theta]", "]"}]}], "+", 
      RowBox[{"\[Epsilon]", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], ",", 
     RowBox[{"\[Epsilon]", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"\[Omega]", " ", "t"}], "+", "\[Phi]"}], "]"}]}], ",", 
     RowBox[{"h1", " ", 
      RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"H2", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"h2", " ", 
       RowBox[{"Sin", "[", "\[Theta]", "]"}]}], "+", 
      RowBox[{"\[Epsilon]", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], ",", 
     RowBox[{"\[Epsilon]", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"\[Omega]", " ", "t"}], "+", "\[Phi]"}], "]"}]}], ",", 
     RowBox[{"h2", " ", 
      RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"H3", "[", "t_", "]"}], ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"h3", " ", 
       RowBox[{"Sin", "[", "\[Theta]", "]"}]}], "+", 
      RowBox[{"\[Epsilon]", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], ",", 
     RowBox[{"\[Epsilon]", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"\[Omega]", " ", "t"}], "+", "\[Phi]"}], "]"}]}], ",", 
     RowBox[{"h3", " ", 
      RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LLG1", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mA", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H1", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mA", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mB", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mAx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mA", "'"}], "[", "t", "]"}]}]}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H1", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mB", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mA", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mBx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mB", "'"}], "[", "t", "]"}]}]}]}]}], "]"}]}], "}"}], "//",
       "Flatten"}], ")"}], "~", "Join", "~", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"mAx", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mAy", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mAz", "[", "0", "]"}], "\[Equal]", 
       SqrtBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]m", "^", "2"}]}]]}], ",", 
      RowBox[{
       RowBox[{"mBx", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mBy", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mBz", "[", "0", "]"}], "\[Equal]", 
       RowBox[{"-", 
        SqrtBox[
         RowBox[{"1", "-", 
          RowBox[{"\[Delta]m", "^", "2"}]}]]}]}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LLG2", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mA", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H2", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mA", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mB", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mAx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mA", "'"}], "[", "t", "]"}]}]}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H2", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mB", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mA", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mBx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mB", "'"}], "[", "t", "]"}]}]}]}]}], "]"}]}], "}"}], "//",
       "Flatten"}], ")"}], "~", "Join", "~", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"mAx", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mAy", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mAz", "[", "0", "]"}], "\[Equal]", 
       SqrtBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]m", "^", "2"}]}]]}], ",", 
      RowBox[{
       RowBox[{"mBx", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mBy", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mBz", "[", "0", "]"}], "\[Equal]", 
       RowBox[{"-", 
        SqrtBox[
         RowBox[{"1", "-", 
          RowBox[{"\[Delta]m", "^", "2"}]}]]}]}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LLG3", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mA", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H3", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mA", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mB", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mAz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mAx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mA", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mA", "'"}], "[", "t", "]"}]}]}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Thread", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "'"}], "[", "t", "]"}], "\[Equal]", 
          RowBox[{
           RowBox[{
            RowBox[{"H3", "[", "t", "]"}], "\[Cross]", 
            RowBox[{"mB", "[", "t", "]"}]}], "+", 
           RowBox[{"\[Omega]E", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"mA", "[", "t", "]"}]}]}], "+", 
           RowBox[{"\[Omega]K", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", 
               RowBox[{"mBz", "[", "t", "]"}]}], "}"}]}]}], "+", 
           RowBox[{"\[Omega]a", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"mBx", "[", "t", "]"}], ",", "0", ",", "0"}], 
              "}"}]}]}], "+", 
           RowBox[{"\[Alpha]", " ", 
            RowBox[{
             RowBox[{"mB", "[", "t", "]"}], "\[Cross]", 
             RowBox[{
              RowBox[{"mB", "'"}], "[", "t", "]"}]}]}]}]}], "]"}]}], "}"}], "//",
       "Flatten"}], ")"}], "~", "Join", "~", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"mAx", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mAy", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mAz", "[", "0", "]"}], "\[Equal]", 
       SqrtBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]m", "^", "2"}]}]]}], ",", 
      RowBox[{
       RowBox[{"mBx", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"mBy", "[", "0", "]"}], "\[Equal]", "\[Delta]m"}], ",", 
      RowBox[{
       RowBox[{"mBz", "[", "0", "]"}], "\[Equal]", 
       RowBox[{"-", 
        SqrtBox[
         RowBox[{"1", "-", 
          RowBox[{"\[Delta]m", "^", "2"}]}]]}]}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LFM", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"\[Phi]", "=", 
       RowBox[{"2", "\[Pi]", "*", "0.01", "*", "j"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"sol", "=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{"LLG1", ",", "variables", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"MaxStepSize", "\[Rule]", "0.001"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Is", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "2", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "1", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}]}]}], ")"}], 
        "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "5", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "4", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "4", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "5", "]"}], "]"}], ",", "t"}], "]"}]}]}], 
         ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{"Is", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"tmax", "/", "2"}], ",", "tmax"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "100"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"SF", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"\[Phi]", "=", 
       RowBox[{"2", "\[Pi]", "*", "0.01", "*", "j"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"sol", "=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{"LLG2", ",", "variables", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"MaxStepSize", "\[Rule]", "0.001"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Is", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "2", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "1", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}]}]}], ")"}], 
        "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "5", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "4", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "4", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "5", "]"}], "]"}], ",", "t"}], "]"}]}]}], 
         ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{"Is", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"tmax", "/", "2"}], ",", "tmax"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "100"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"QFM", "=", 
   RowBox[{"ParallelTable", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"\[Phi]", "=", 
       RowBox[{"2", "\[Pi]", "*", "0.01", "*", "j"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"sol", "=", 
       RowBox[{"NDSolveValue", "[", 
        RowBox[{"LLG3", ",", "variables", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"MaxStepSize", "\[Rule]", "0.001"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Is", "=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "2", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "1", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "1", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "2", "]"}], "]"}], ",", "t"}], "]"}]}]}], ")"}], 
        "+", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "5", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "4", "]"}], "]"}], ",", "t"}], "]"}]}], "-", 
          RowBox[{
           RowBox[{"sol", "[", 
            RowBox[{"[", "4", "]"}], "]"}], 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"sol", "[", 
              RowBox[{"[", "5", "]"}], "]"}], ",", "t"}], "]"}]}]}], 
         ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"NIntegrate", "[", 
       RowBox[{"Is", ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", 
          RowBox[{"tmax", "/", "2"}], ",", "tmax"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "100"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"LFM", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Black", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"SF", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Magenta", ",", "Dashed"}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", 
  RowBox[{"QFM", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Darker", "[", "Cyan", "]"}], ",", "Thick", ",", "DotDashed"}], 
     "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.784667182320609*^9, 3.784667253820386*^9}, {
   3.7846752636144543`*^9, 3.7846752911270065`*^9}, {3.7846754138453445`*^9, 
   3.784675416554332*^9}, {3.7846754513135805`*^9, 3.784675487638877*^9}, {
   3.784675525588106*^9, 3.7846755320339475`*^9}, {3.7846756195391903`*^9, 
   3.78467576093281*^9}, {3.7846757941111374`*^9, 3.7846758357232466`*^9}, 
   3.7846922901004915`*^9, {3.7847064725978093`*^9, 3.7847065950360193`*^9}, {
   3.7847747224921074`*^9, 3.784774749186576*^9}, {3.784774789479787*^9, 
   3.784774823229568*^9}, {3.7847773896573825`*^9, 3.7847773994838457`*^9}, {
   3.7847805174268575`*^9, 3.7847805203331146`*^9}, {3.7848520599651976`*^9, 
   3.784852066561028*^9}},
 CellLabel->"In[18]:=",ExpressionUUID->"a354ad9c-c7e5-4ac0-9224-66c264ac14ac"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {GrayLevel[0], PointSize[0.009166666666666668], Thickness[Large], 
     LineBox[CompressedData["
1:eJw11QlMk3cYx/EKOM4wIAyRo3IJ0gFyidz9UShXOdrSQjm2sSmwTblkMkFC
qmI4JjBgDpidUxcqcxAQJolLlDoy7iGSsXlxhRmVw7kxFZDhXPq8TZrm02+f
5M2b9//U9oNccYYGi8XKeP3+/1P9esrV54V7vxEboaIv8NkT536reMY6sE3b
UKTzGRvBcrn21Z3djE2R2WpyRmHI2BznXPU+HZgLJ1shpyxvruw7xjaYOqmM
0c1kbIetn4y4N21j7IDNQ9WjRdf5ZEecrFx4ppXKeBfWM1PGm5fDyBzMyhJL
LxQydoHC7+5Q/fNQshvWJ5xvPMxm7I6D7RKB0xSP7AGlwZeXZ8IYe8Ky0eur
sdYQshfmJW32lizG3ljKNNSdlYC8B2nRxUpEcck+2By+amE+G0Tei/tJLeLO
3ECyL06ZddQcXfMn+6HJZlpjpdSP7I8w9k8pp7f4kgOQYvGr14PjPuRA9CbX
1mZp7SEHYWS/jHuh3IscjJW8fQILQ08yF/PsYZm42Z0MtI5ratVL3dSWA5zb
bl3L5W9TD8Eqfx/77Ogu6iGocZ49fsjRkToPVfOZHOcme+o8dNd2qrIcbKmH
omer8nDBGWvqoehIv/LSzHU79TAslUa9n6TxFvUwvOdREle5w5g6H6/cAlS5
xvrU+TDOf7cn+IQW9XC0uF65963Pv73qHg7PB0veit6/1WZF4MhHVburZfPU
I1DuVrSc7a6iHgn22Ngf1W3jXHWPxHNRcdnkjTm1WVGw9067tODyiHoULAYi
rfXqFqlHo/lmTNbDBbI8GhncmN4P/R9TFyCDbemhPD5PXQDrzbtbhhsmqcfg
nfP2PxuNKej6YsA1Xv2tXXOGri8Wt1mP8tNLnlKPxa3UjiWJYp16HPSud4v1
9TXp/sQhcq2neyBHl+5PPFbEQb4SPyPq8Yh1agpa1DalLkQb9/IP2WxztSFE
mYlsNbnIin4vhLvOYPN6pY3aKiH6/Pffu/SPHc2LMG23IRUW76R5EbJqTuUL
zJnnRYTGSu2DiyMcmhehXUfb70C9K82LUVjkEFjsTc8jxFjuPhsgUnrQvBhm
B4YUPGt6nlWvfz9Vypto9Kb5BEy0VBQeM6XzgAS8yfXKGa/bS/MJGOzKHeEa
0nlSJWB0jb1dXMGcNwnO96/wqjcDaF4CQdznm0dz6bzKJQhp+yWw404wzUvw
+7n712asmfMuRX9qn6HTMBlS4PCL+tU82g9yKaS3rhmMG9M+UUmhFIf9+eR7
Zr8k4ouX/IbBYNpHSITLtsqF0SGyPBFzH3MGCmJon6kSYeLgIQzqZ/ZbElgh
Fwu6fGgfIgmm06FxGV+T5UmIPPai6q81sioJFZoGP24IwmleBv3Jum98GsiQ
4WLC40XcJMtl2MF/5s9h0T5XyXC6sLPvyE5mvyejvWRd24ZLRjKuRqj/P/4D
L+7a2Q==
      "]]}}, {}, {}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 101.}, {-46.61806561515285, 0.8724247949198684}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.784667546363476*^9, 3.7846755259451294`*^9, 
  3.784676410627068*^9, 3.7846928767140584`*^9, 3.7847074775467396`*^9, 
  3.7847774968980513`*^9, 3.784781479060573*^9, 3.7848530057913027`*^9},
 CellLabel->"Out[31]=",ExpressionUUID->"ed230f57-7ca6-41e9-b802-7ef910f51289"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[1, 0, 1], PointSize[0.009166666666666668], AbsoluteThickness[
     1.6], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJw1lXtMU2cYh0s7BMlEqN24FBwUxEnAWUCEMewPsNyvpbQF8cIAgTJlmYTb
tlp1DpmBcVFAwBDiojGEKX9gkYnUcBk6w1h3EReGCA6HqCAYGOOyG997kubk
6XOec75+eU/q+H6OLJ3L4XDS//38d147ZiRH9+rH7i0J9ewLKPfd2l4+SmyK
i8f42YJeYgsYq2IKeM3EApTNF5ZPfEFsjfz6xC02+4ntcEfdPRTvS+yANK6l
wXiSWASl7Y3Urut2jJ3RFDM+NWVD7ILW9Qs5Xjpbxm/jtCGys+E+sStq1SNH
FAt0PzeYBGae98mgfjt+HPjj3OM58juwYvJlqX81sRi9JQvPxz8i9kDOW5qi
po+JPSEXbLQzXCD2gteNq22hA8Q7UeUyb8fl0vO84eR8byVnB/EucNWL22RK
Yh8YS0K9e/OJfSE+VWP1oor4XXSP/5Id2ULshxEIe9xvEr+HC5UD0/wsYn+E
j9wVen9C69mNX8XfJz6Lo/2RYNtx/VzRkjVj4GXZFcveXKs11gL1Ys1+bdsb
zAdAfVLhJGnfxHwANKqJJ13FlswHon/ZrPl6rTnzgciymE9qzDdjPgjLnLzZ
DZPrmA9Cs8rrIGZ4zO8B/1xui1uBEfN7IB/7OnnCsNq15qWQmvqN/pyytMZa
KXjKDUKp9k/mg3EwRayetJ5nPhjaz6ce2tm9Yj4EHsKlS8maWeZD4H45ieuz
6yXzoXjs9+ncWcwwHwr/rSJemts082EY233X77L6BfNh4Dnm1zn5MuaE49hw
wO+aGvLhqEttDRsspj4CNZyN6rRiun8EMvpWl9fV0fMjEVwdxcsymmM+EjU6
s8MHemn9UZh3FVZemaHfF4XOlMbR6YpF5qPxwOREgFvnMvPR0B8KTu8K5rD9
jUG0bVtfXhWX7W8MPI9/VdWbZMx8LHCHy+u4ZrrGiIVOk/0Xv+p1dn0sjnTP
5RqnW6yxPha31D19Twv4rI9DS7tRYN4BAevjsPVhRnqJ+Zusj0P6+Q8ixGfY
fOnj4FLswzUfofmTwd7K6vk35mw+IcNQU43buIjNr1aGhrMm/PWb2XzrZdDx
hwVJHHvWx2P08NDgt7cZIx72Pz3adMJ+M+vjcanUJpBnylgfj5KTv30XM0q9
HK9KK7p17dTLYfbEpfBZBWOtHPkmhsSsbMZ6OdxFP1T6SKlPwGfTZldfc6A+
Af1yo1DRCluvNgGaB82pniO0/gQ0FE/uO9VD76sCtU/vZ/dfYwwF+hv6LV0u
Uq/AcOY7HlvqqVdA6R60c6aaeiV0NsnlYWeoV6KMZ5EgOUq9EqLC0VbOIeqV
6LDNuy1JpV4F59zGm0Vx1Kvw94d7F4+FUK+C1OC4+siPehWEsytWUlfqE+Ew
2N96WkB9InwyO/7///gHwDbYMA==
      "]]}}, {}, {}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 101.}, {-7.528816519411547, 0}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.784667546363476*^9, 3.7846755259451294`*^9, 
  3.784676410627068*^9, 3.7846928767140584`*^9, 3.7847074775467396`*^9, 
  3.7847774968980513`*^9, 3.784781479060573*^9, 3.7848530058694315`*^9},
 CellLabel->"Out[32]=",ExpressionUUID->"e24ac485-f45d-4cd6-97ba-e431b26a6a2b"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[0, 
      NCache[
       Rational[2, 3], 0.6666666666666666], 
      NCache[
       Rational[2, 3], 0.6666666666666666]], PointSize[0.009166666666666668], 
     Thickness[Large], Dashing[{0, Small, Small, Small}], 
     LineBox[CompressedData["
1:eJw11XlQlHUcx/EVImigBbkl5FA0g1hMCMFF+IDcy7EHLCwKrCTZkqQTGQge
D0eoeEwgKWjpwJCT5WCAxAgjPhxqAxSIIS0QyCTKGagQJofV7PfZmZ2d176f
549n5/f9rn3SXmmyFo/HS/7v/f+n5jXjMygfUy2F2bP0BZgC76fCUM56+Kiv
vdLZh7MRyovWfPnQibMpmvIT1TFGnC3RFfJysGDSjmyNY/tVN8+ynO3gkfXr
ztyTnNcgt71+1W4xZwd8EjyuiDfgvB6HtAWt/Gxb8gZUj6fH2fbZkB2R6MYf
lTpwfhdxG+aHrJWryQL4H45QGp+2Jm9EYVz9ur7Kt8jvYffNAyqPFivyJkha
7yQZdqwiu+JU4MiCss2S7Abn/hKZ+JYF+X18p6hpsL9qTnZHVq/LVc8iM/Jm
PMhb52WVZkr2QEpJWV+C2ITsCUM2YvtPTsbkLXBzupPZ//pKshDXuyf3WN41
JHthR/jUK97RN8lbwS40PGoS65O9ob/2RerBt98g+6DjmuteE74uGThpqF5v
p6ujMQOEtvUKJky1qftC6Guep/JcQd0X572UoszJ5Vua7oeocdH5tNIFjRk/
2Ix/43Y/6wX1bdC6MKXl/NUc9W1IMpX6VDx+Rt0fGY4nP1zaN0PdH/NqQU6C
apJ6AHbt35Uw8sEo9QDMBNVeKbv3iHogiq1GxhOqh6kHYjGwezhVf5B6EDo2
Tv8YeFZNPQg+Y5u+tnjtAfVgOC6faC6X3acejGnLiYzK3HvUQ2AyfNra+EwX
9RB8u6fil0aGzAtFhll7WquIrmdCkd8Jsex5N3URVHzGu/7zHuoilBe2lu3s
+Z16GJxdFh428/+gHoZmudBKsJKehxcOp8QCYZcuPS8TDqWLs+hpyhPqEWjW
abez2TJBPQIpF4ob6z6bph6J08LHdTMG9HszkWAuS23di2epizG6Qns2MXJe
Y4jReDGhX7n1JV0vBm8wYHFgx5LGrBj6UduD+HIenQ8Jckr8s46YamkMCWIy
La6kj9L5YSQYMG9pSu2i88VKkN3dm3yujTt/UiT8PHFo9Dc6n5Cir/DPHrMp
Or+MFF6u6poSIz7dL8WrBoFBvI0R3S/Dqc4yo/DNNB+Q4a+kUctLEpofRobV
6UNOnR/TfLEyiLSfXFLkcfMXBf0CSVtTKc0nonDmtt3zc9/T/DJROLYsiD5c
S/PNRiG/Qp3jeoOb/2gUxZofr7lO+wHRqFOluteupH3CRGOfbvCR0lNkNhrs
pwFuR1Xc/pFjwPLGcaUn7SfIUSYRFCoWyYwcJ+5mS8eqaJ+xckiLu+KzFdz+
i8HtloprPzwjIwYujS4O2d60L5kYvCPsrTY/SGZjYCKd+3umitunsQirUevN
DJERi6pCP71/dGh/M7HI7/a6OLeWzMZiVivwQLcnt98V0L6s8vwiiAwF6sM0
/x//Am8Z5yU=
      "]]}}, {}, {}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 101.}, {-8.65037089322453, 0}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.784667546363476*^9, 3.7846755259451294`*^9, 
  3.784676410627068*^9, 3.7846928767140584`*^9, 3.7847074775467396`*^9, 
  3.7847774968980513`*^9, 3.784781479060573*^9, 3.7848530059475574`*^9},
 CellLabel->"Out[33]=",ExpressionUUID->"1877d71d-6dff-4464-af06-c8584347e07c"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"LFM", ",", "SF", ",", "QFM"}], "}"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"FrameStyle", "\[Rule]", "Medium"}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"Black", ",", "Thick"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Magenta", ",", "Dashed"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Darker", "[", "Cyan", "]"}], ",", "Thick", ",", 
        "DotDashed"}], "}"}]}], "}"}]}]}], "]"}]], "Input",
 CellLabel->"In[34]:=",ExpressionUUID->"9355ab40-618f-4804-af9b-e3c06963aebd"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {GrayLevel[0], PointSize[0.007333333333333334], Thickness[Large], 
     LineBox[CompressedData["
1:eJw11QlMk3cYx/EKOM4wIAyRo3IJ0gFyidz9UShXOdrSQjm2sSmwTblkMkFC
qmI4JjBgDpidUxcqcxAQJolLlDoy7iGSsXlxhRmVw7kxFZDhXPq8TZrm02+f
5M2b9//U9oNccYYGi8XKeP3+/1P9esrV54V7vxEboaIv8NkT536reMY6sE3b
UKTzGRvBcrn21Z3djE2R2WpyRmHI2BznXPU+HZgLJ1shpyxvruw7xjaYOqmM
0c1kbIetn4y4N21j7IDNQ9WjRdf5ZEecrFx4ppXKeBfWM1PGm5fDyBzMyhJL
LxQydoHC7+5Q/fNQshvWJ5xvPMxm7I6D7RKB0xSP7AGlwZeXZ8IYe8Ky0eur
sdYQshfmJW32lizG3ljKNNSdlYC8B2nRxUpEcck+2By+amE+G0Tei/tJLeLO
3ECyL06ZddQcXfMn+6HJZlpjpdSP7I8w9k8pp7f4kgOQYvGr14PjPuRA9CbX
1mZp7SEHYWS/jHuh3IscjJW8fQILQ08yF/PsYZm42Z0MtI5ratVL3dSWA5zb
bl3L5W9TD8Eqfx/77Ogu6iGocZ49fsjRkToPVfOZHOcme+o8dNd2qrIcbKmH
omer8nDBGWvqoehIv/LSzHU79TAslUa9n6TxFvUwvOdREle5w5g6H6/cAlS5
xvrU+TDOf7cn+IQW9XC0uF65963Pv73qHg7PB0veit6/1WZF4MhHVburZfPU
I1DuVrSc7a6iHgn22Ngf1W3jXHWPxHNRcdnkjTm1WVGw9067tODyiHoULAYi
rfXqFqlHo/lmTNbDBbI8GhncmN4P/R9TFyCDbemhPD5PXQDrzbtbhhsmqcfg
nfP2PxuNKej6YsA1Xv2tXXOGri8Wt1mP8tNLnlKPxa3UjiWJYp16HPSud4v1
9TXp/sQhcq2neyBHl+5PPFbEQb4SPyPq8Yh1agpa1DalLkQb9/IP2WxztSFE
mYlsNbnIin4vhLvOYPN6pY3aKiH6/Pffu/SPHc2LMG23IRUW76R5EbJqTuUL
zJnnRYTGSu2DiyMcmhehXUfb70C9K82LUVjkEFjsTc8jxFjuPhsgUnrQvBhm
B4YUPGt6nlWvfz9Vypto9Kb5BEy0VBQeM6XzgAS8yfXKGa/bS/MJGOzKHeEa
0nlSJWB0jb1dXMGcNwnO96/wqjcDaF4CQdznm0dz6bzKJQhp+yWw404wzUvw
+7n712asmfMuRX9qn6HTMBlS4PCL+tU82g9yKaS3rhmMG9M+UUmhFIf9+eR7
Zr8k4ouX/IbBYNpHSITLtsqF0SGyPBFzH3MGCmJon6kSYeLgIQzqZ/ZbElgh
Fwu6fGgfIgmm06FxGV+T5UmIPPai6q81sioJFZoGP24IwmleBv3Jum98GsiQ
4WLC40XcJMtl2MF/5s9h0T5XyXC6sLPvyE5mvyejvWRd24ZLRjKuRqj/P/4D
L+7a2Q==
      "]]}, 
    {RGBColor[1, 0, 1], PointSize[0.007333333333333334], AbsoluteThickness[
     1.6], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJw1lXtMU2cYh0s7BMlEqN24FBwUxEnAWUCEMewPsNyvpbQF8cIAgTJlmYTb
tlp1DpmBcVFAwBDiojGEKX9gkYnUcBk6w1h3EReGCA6HqCAYGOOyG997kubk
6XOec75+eU/q+H6OLJ3L4XDS//38d147ZiRH9+rH7i0J9ewLKPfd2l4+SmyK
i8f42YJeYgsYq2IKeM3EApTNF5ZPfEFsjfz6xC02+4ntcEfdPRTvS+yANK6l
wXiSWASl7Y3Urut2jJ3RFDM+NWVD7ILW9Qs5Xjpbxm/jtCGys+E+sStq1SNH
FAt0PzeYBGae98mgfjt+HPjj3OM58juwYvJlqX81sRi9JQvPxz8i9kDOW5qi
po+JPSEXbLQzXCD2gteNq22hA8Q7UeUyb8fl0vO84eR8byVnB/EucNWL22RK
Yh8YS0K9e/OJfSE+VWP1oor4XXSP/5Id2ULshxEIe9xvEr+HC5UD0/wsYn+E
j9wVen9C69mNX8XfJz6Lo/2RYNtx/VzRkjVj4GXZFcveXKs11gL1Ys1+bdsb
zAdAfVLhJGnfxHwANKqJJ13FlswHon/ZrPl6rTnzgciymE9qzDdjPgjLnLzZ
DZPrmA9Cs8rrIGZ4zO8B/1xui1uBEfN7IB/7OnnCsNq15qWQmvqN/pyytMZa
KXjKDUKp9k/mg3EwRayetJ5nPhjaz6ce2tm9Yj4EHsKlS8maWeZD4H45ieuz
6yXzoXjs9+ncWcwwHwr/rSJemts082EY233X77L6BfNh4Dnm1zn5MuaE49hw
wO+aGvLhqEttDRsspj4CNZyN6rRiun8EMvpWl9fV0fMjEVwdxcsymmM+EjU6
s8MHemn9UZh3FVZemaHfF4XOlMbR6YpF5qPxwOREgFvnMvPR0B8KTu8K5rD9
jUG0bVtfXhWX7W8MPI9/VdWbZMx8LHCHy+u4ZrrGiIVOk/0Xv+p1dn0sjnTP
5RqnW6yxPha31D19Twv4rI9DS7tRYN4BAevjsPVhRnqJ+Zusj0P6+Q8ixGfY
fOnj4FLswzUfofmTwd7K6vk35mw+IcNQU43buIjNr1aGhrMm/PWb2XzrZdDx
hwVJHHvWx2P08NDgt7cZIx72Pz3adMJ+M+vjcanUJpBnylgfj5KTv30XM0q9
HK9KK7p17dTLYfbEpfBZBWOtHPkmhsSsbMZ6OdxFP1T6SKlPwGfTZldfc6A+
Af1yo1DRCluvNgGaB82pniO0/gQ0FE/uO9VD76sCtU/vZ/dfYwwF+hv6LV0u
Uq/AcOY7HlvqqVdA6R60c6aaeiV0NsnlYWeoV6KMZ5EgOUq9EqLC0VbOIeqV
6LDNuy1JpV4F59zGm0Vx1Kvw94d7F4+FUK+C1OC4+siPehWEsytWUlfqE+Ew
2N96WkB9InwyO/7///gHwDbYMA==
      "]]}, 
    {RGBColor[0, 
      NCache[
       Rational[2, 3], 0.6666666666666666], 
      NCache[
       Rational[2, 3], 0.6666666666666666]], PointSize[0.007333333333333334], 
     Thickness[Large], Dashing[{0, Small, Small, Small}], 
     LineBox[CompressedData["
1:eJw11XlQlHUcx/EVImigBbkl5FA0g1hMCMFF+IDcy7EHLCwKrCTZkqQTGQge
D0eoeEwgKWjpwJCT5WCAxAgjPhxqAxSIIS0QyCTKGagQJofV7PfZmZ2d176f
549n5/f9rn3SXmmyFo/HS/7v/f+n5jXjMygfUy2F2bP0BZgC76fCUM56+Kiv
vdLZh7MRyovWfPnQibMpmvIT1TFGnC3RFfJysGDSjmyNY/tVN8+ynO3gkfXr
ztyTnNcgt71+1W4xZwd8EjyuiDfgvB6HtAWt/Gxb8gZUj6fH2fbZkB2R6MYf
lTpwfhdxG+aHrJWryQL4H45QGp+2Jm9EYVz9ur7Kt8jvYffNAyqPFivyJkha
7yQZdqwiu+JU4MiCss2S7Abn/hKZ+JYF+X18p6hpsL9qTnZHVq/LVc8iM/Jm
PMhb52WVZkr2QEpJWV+C2ITsCUM2YvtPTsbkLXBzupPZ//pKshDXuyf3WN41
JHthR/jUK97RN8lbwS40PGoS65O9ob/2RerBt98g+6DjmuteE74uGThpqF5v
p6ujMQOEtvUKJky1qftC6Guep/JcQd0X572UoszJ5Vua7oeocdH5tNIFjRk/
2Ix/43Y/6wX1bdC6MKXl/NUc9W1IMpX6VDx+Rt0fGY4nP1zaN0PdH/NqQU6C
apJ6AHbt35Uw8sEo9QDMBNVeKbv3iHogiq1GxhOqh6kHYjGwezhVf5B6EDo2
Tv8YeFZNPQg+Y5u+tnjtAfVgOC6faC6X3acejGnLiYzK3HvUQ2AyfNra+EwX
9RB8u6fil0aGzAtFhll7WquIrmdCkd8Jsex5N3URVHzGu/7zHuoilBe2lu3s
+Z16GJxdFh428/+gHoZmudBKsJKehxcOp8QCYZcuPS8TDqWLs+hpyhPqEWjW
abez2TJBPQIpF4ob6z6bph6J08LHdTMG9HszkWAuS23di2epizG6Qns2MXJe
Y4jReDGhX7n1JV0vBm8wYHFgx5LGrBj6UduD+HIenQ8Jckr8s46YamkMCWIy
La6kj9L5YSQYMG9pSu2i88VKkN3dm3yujTt/UiT8PHFo9Dc6n5Cir/DPHrMp
Or+MFF6u6poSIz7dL8WrBoFBvI0R3S/Dqc4yo/DNNB+Q4a+kUctLEpofRobV
6UNOnR/TfLEyiLSfXFLkcfMXBf0CSVtTKc0nonDmtt3zc9/T/DJROLYsiD5c
S/PNRiG/Qp3jeoOb/2gUxZofr7lO+wHRqFOluteupH3CRGOfbvCR0lNkNhrs
pwFuR1Xc/pFjwPLGcaUn7SfIUSYRFCoWyYwcJ+5mS8eqaJ+xckiLu+KzFdz+
i8HtloprPzwjIwYujS4O2d60L5kYvCPsrTY/SGZjYCKd+3umitunsQirUevN
DJERi6pCP71/dGh/M7HI7/a6OLeWzMZiVivwQLcnt98V0L6s8vwiiAwF6sM0
/x//Am8Z5yU=
      "]]}}, {}, {}, {}, {}},
  AspectRatio->1,
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0., 0},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameStyle->Medium,
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0., 101.}, {-46.61806561515285, 0.8724247949198684}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.7847815765408525`*^9, 3.7848533145366817`*^9},
 CellLabel->"Out[34]=",ExpressionUUID->"83f4244b-7cae-4149-a7d7-4bb5dd536b2b"]
}, Open  ]]
},
WindowSize->{1206, 1021},
WindowMargins->{{Automatic, 34}, {Automatic, 0}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 14310, 342, 500, "Input",ExpressionUUID->"5981de4b-bc46-466a-bea8-49eb3958bb2a"],
Cell[14893, 366, 15265, 266, 238, "Output",ExpressionUUID->"82d99a53-3f41-4f55-9d31-2bd4476e5722"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30195, 637, 1042, 17, 73, "Input",ExpressionUUID->"b6431c20-3ae4-4860-ab7e-b875fe6ca575"],
Cell[31240, 656, 726, 12, 34, "Output",ExpressionUUID->"ed1c2fcc-0ad6-41e5-83e7-f66d6fb17979"],
Cell[31969, 670, 725, 12, 34, "Output",ExpressionUUID->"036d09fd-b7ef-4dd2-9654-e832b4773e44"],
Cell[32697, 684, 724, 12, 34, "Output",ExpressionUUID->"cb8b5f97-1376-45f5-ab4a-9d35c5b21037"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33458, 701, 20132, 558, 682, "Input",ExpressionUUID->"a354ad9c-c7e5-4ac0-9224-66c264ac14ac"],
Cell[53593, 1261, 2887, 62, 244, "Output",ExpressionUUID->"ed230f57-7ca6-41e9-b802-7ef910f51289"],
Cell[56483, 1325, 2865, 61, 249, "Output",ExpressionUUID->"e24ac485-f45d-4cd6-97ba-e431b26a6a2b"],
Cell[59351, 1388, 3031, 67, 249, "Output",ExpressionUUID->"1877d71d-6dff-4464-af06-c8584347e07c"]
}, Open  ]],
Cell[CellGroupData[{
Cell[62419, 1460, 782, 20, 52, "Input",ExpressionUUID->"9355ab40-618f-4804-af9b-e3c06963aebd"],
Cell[63204, 1482, 6059, 121, 370, "Output",ExpressionUUID->"83f4244b-7cae-4149-a7d7-4bb5dd536b2b"]
}, Open  ]]
}
]
*)

